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Abstract 

The development of quantum cosmology, in which Stephen Hawk- 
ing played a crucial role, has frequently encountered substantial con- 
ceptual and technical difficulties related to the problem of time in 
quantum gravity and to general issues concerning the foundations of 
quantum theory. In this contribution to Stephen's 60th Birthday Con- 
ference, I describe some recent work in which the decoherent histories 
approach to quantum theory is used to quantize simple cosmological 
models and perhaps shed some light on some of these difficulties. 
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1 Introduction 

It is a truly great pleasure to be part of Stephen Hawking's 60th birthday 
celebrations. It is now almost twenty years since I commenced my PhD 
studies with Stephen, and I recall that exciting time very fondly. Numerous 
researchers were just beginning to address the conceptual and technical dif- 
ficulties of quantum cosmology, and the confusion and controversy at that 
time led to some very lively debates in which Stephen played a very central 
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role. The atmosphere of this conference has reminded me very pleasantly of 
those early days. 

The work described here concerns some generally interesting issues in 
quantum cosmology. This is the area I started working on with Stephen, and 
its intriguing problems and difficulties have been the inspiration for much 
else I have pursued in physics. Quantum cosmology was born in the 1960's 
with the work of Bryce DeWitt, Charles Misner, John Wheeler and others 
[|, ^7], [3^], but it really took off in the early 1980's with Hartle and Hawking's 



seminal paper, The Wave Function of the Universe fll8| . This development 



was very pertinent since the appearance of the inflationary universe scenario 
at about that time underscored the urgency to acquire an understanding of 
cosmological initial conditions. The new ingredients of quantum cosmology, 
compared to its older 1960s version, were two-fold. First, there was the use of 
Euclidean path integral methods which had grown out of their very successful 
application to black holes physics. Secondly, and perhaps more importantly, 
there was the Hartle-Hawking "no-boundary" proposal. This gave, perhaps 
for the first time, a genuinely quantum gravitational proposal which implied 
cosmological initial conditions. In this proposal, wave functions satisfying 
the Wheeler-DeWitt equation 

H*[h ij \=Q (1) 
are given by a sum over histories expression of the form 

= / P^exp(-J[^]) (2) 
J c 

where the sum is over closed four-geometries whose only boundary is the 
three-surface on which the three-metric hij is specified. (Note also that the 
functional integral is taken to be over a complex contour C 

There were two questions of particular interest in these early days. Firstly, 
what sort of picture can quantum gravity give about the initial singularity? 
And secondly, what initial conditions does the theory imply for classical 
cosmological solutions? In simple, somewhat heuristic models, the answer 
to the first question is that the initial singularity is replaced by some kind 
of quantum tunneling situation. That is, the classical singularity can (at 
least in some models) become a classically forbidden region. The second 
question can be answered by showing how quantum cosmology produces, in 
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a quasiclassical limit, an approximate probability distribution on the space 
of classical cosmological solutions. 

The most pressing questions in those days were to do with cosmological 
issues: Does the no-boundary proposal predict initial conditions which pro- 



duce inflationary solutions and subsequent structure formation |L9|, |I3| ? The 
interpretation associated with these models was heuristic |24], [8|, 10 , . For 
example, one considered the Klein-Gordon flux 

J = i (\I>*V^/ - \I/V\[>*) (3) 

associated with the wave functions of the WKB type, 

^ = Ce iS (4) 

These heuristic methods appeared to be adequate for the questions consid- 
ered at the time. 

Now, however, in 2002, far more is known about the foundations and 
interpretation of quantum theory. Moreover, there has also been a consid- 
erable amount of activity in the experimental tests of foundational ideas. 
As a result of this, the foundations of quantum theory is now more in the 
mainstream of physics, when previously it was regarded as the domain of the 
philosophers. Perhaps as a result of this, there has been a secondary wave of 
much slower activity in quantum cosmology, that aims to look more deeply 
into a whole variety of mathematical and conceptual issues related to the 
application of quantum theory to the universe as a whole. In particular, it 
is now reasonable to ask, how can the heuristic ideas used earlier be derived 
or reconciled with a properly defined interpretational structure for quantum 
theory of models described by a Wheeler-DeWitt equation? For example, in 
standard quantum theory, probabilities generally have the form 

P = Tr (Pp) (5) 

where p is the density operator, P is a projection operator and the trace 
is over a complete set of states (and so assumes a Hilbert space of states 
with a suitably defined inner product). How can the predictions of quantum 
cosmology be reconciled with a formula of this type? 



Related to these issues is the notorious problem of time [24, 21, H[. The 



Wheeler-DeWitt equation, in simple minisuperspace models, has the form, 

= (-V 2 + U)V = (6) 
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where V 2 is a d'Alembertian type operator. The equation is of the form of 
a Klein-Gordon equation in a general curved spacetime background with a 
space-time dependent mass term. There is no single variable to play the role 
of time, nor is there the possibility of splitting the solutions into positive 
and negative frequency. In the usual Klein-Gordon equation one then resorts 
to second quantization, but the Wheeler-DeWitt equation is in some sense 
already the second quantized theory, so one has to face the issue of the lack 
of a time coordinate more squarely. 

In this contribution to the conference, I will describe how the decoherent 
histories approach may be used to provide a quantization of simple minisu- 
perspace models, perhaps avoiding some of the serious difficulties outlined 
above, and agreeing with the heuristic methods used earlier on. 

I begin by making two simple observations about the Wheeler-DeWitt 
equation that are generally relevant to the discussion. The first is that, at 
least classical speaking, the constraint equation, 

H = f ab PaPb + V = (7) 

is related to reparametrization invariance. (This is the leftover of four- 
dimensional diffeomorphism invariance after the restriction to minisuper- 
space). In keeping with Dirac's general ideas about quantizing constrained 
systems, we therefore look for observables- quantities which commute with 
the constraint. There has been a certain amount of debate about this issue 



in the context of general relativity, but this will not affect us here [y, [33[ . 

Second, as observed by Barbour [[1], |2], by analogy with Mott's 1929 cal- 
culation of tracks in a cloud chamber p9| , there is a natural association 
between the Wheeler-DeWitt equation and the emergence of classical trajec- 
tories. Mott asked why the outgoing spherical wave associated with alpha 
decay produced a straight line track in a cloud chamber. By considering the 
Hamiltonian of the decaying atom interacting with ionizing particles in the 
cloud chamber he showed that the ionized particles with high probability lie 
along a straight line. Interestingly, and as Barbour observed, Mott actually 
solved the time-independent Schrodinger equation, 

= (8) 

rather than the time-dependent one, to obtain the tracks. This is therefore 
an interesting analogy for the Wheeler-DeWitt equation where, in the quasi- 
classical limit one expects to obtain emergent classical trajectories. (In fact, 
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a model of the Wheeler-DeWitt equation including explicit detectors, and 
exhibiting these features has been constructed |J. Some related approaches 
to timeless models are Refs. |5T] , |52fl.) 



With these preliminaries out of the way, we now focus on the following 
simple question: Suppose we have an n dimensional configuration space with 
coordinates x = (x\, x 2 , ■ ■ -x n ), and suppose the wave function of the system 
is in an eigenstate of the Hamiltonian, 

HV(xi,x 2 , ■ ■ -Xn) = E^(xi,x 2 , ■ ■ -x n ) (9) 

What is the probability of finding the system in a series of regions of config- 
uration space Ai, A 2 , • ■ ■ A N without reference to time? 

This question will form the main focus of the rest of this paper. We first 
consider the classical case, and then the decoherent histories analysis. 



2 The Classical Case 

We begin by considering the classical case which contains many almost all 



the key features of the problem. (We follow the treatment of Ref.[|15[ quite 
closely). For simplicity we will concentrate on the case of a single region of 
configuration space A. 

We will consider a classical system described by a 2n- dimensional phase 
space, with coordinates and momenta (x, p) = (xk,Pk), and Hamiltonian 

H = m + y M 

More generally, we are interested in a system for which the kinetic part of 
the Hamiltonian has the form g h i(x)pkPj, where (^(x) is an inverse metric 
of hyperbolic signature. Most minisuperspace models in quantum cosmology 
have a Hamiltonian of this form. However, the focus of this paper is the 
timelessness of the system, and the form of the configuration space metric 
turns out to be unimportant. So for simplicity, we will concentrate on the 
form Eq.fllOD. 

We assume that there is a classical phase space distribution function 
w(p, x), which is normalized according to 

J d n p d n x w(p, x) = 1 (11) 
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and obeys the evolution equation 

dw ^ / p k dw dV dw\ 

where { , } denotes the Poisson bracket. The interesting case is that in which 
w is the classical analogue of an energy eigenstate, in which case dw/dt = 0, 
so the evolution equation is simply 

{H,w} = (13) 

It follows that 

W (p c/ (t),x d (t))= W (p(0),x(0)) (14) 

where p d (t), x c/ (t) are the classical solutions with initial data p(0), x(0), so w 
is constant along the classical orbits. (The normalization of w then becomes 
an issue if the classical orbits are infinite, but we will return to this in the 
quantum case discussed below). 

Given a set of classical solutions (p d (t), x d (t)), and a phase space dis- 
tribution function w, we are interested in the probability that a classical 
solution will pass through a region A of configuration space. We construct 
this as follows. First of all we introduce the characteristic function of the 
region A, 

W = {o itwt. < 15 > 

To see whether the classical trajectory x. cl (t) intersects this region, consider 
the phase space function 

A(x, Po ,x ) = r dt 5 w (x-x ci (f)) (16) 

J — oo 

(In the case of periodic classical orbits, the range of t is taken to be equal to 
the period). This function is positive for points x on the classical trajectory 
labeled by Po,x and zero otherwise. Hence intersection of the classical 
trajectory with the region A means, 

J d n x / A (x) y°° dt 5 (n) (x - x d (t)) > (17) 

Or equivalently, that 

r dt / A (x d (t)) > o (is) 

J — oo 
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This quantity is essentially the amount of parameter time the trajectory 
spends in the region A. We may now write down the probability for a 
classical trajectory entering the region A. It is, 

Pa — J d n p d n x w(p , x ) 6 QT dt / A (x d (t)) - e) (19) 

In this construction, e is a small positive number that is eventually sent to 
zero, and is included to avoid possible ambiguities in the ^-function at zero 
argument. The ^-function ensures that the phase space integral is over all 
initial data whose corresponding classical trajectories spend a time greater 
than e in the region A. 

The classical solution x d (i) depends on some fiducial initial coordinates 
and momenta, x and po, say In the case of a free particle, for example, 

x^)=x + ^ (20) 

The construction is independent of the choice of fiducial initial points. If we 
shift x , Po along the classical trajectories, the measure, phase space distri- 
bution function w and the function are all invariant. Hence the integral 
over xo, po is effectively a sum over classical trajectories. The shift along the 
classical trajectories may also be thought of as a reparametrization, and the 
quantity (|19"D is in fact a reparametrization-invariant expression of the notion 
of a classical trajectory. This means that the probability (|T9"D has the form 



of a phase space overlap of the "state" with a reparametrization-invariant 
operator. 

It is useful also to write this result in a different form, which will be 
more relevant to the results we get in the quantum theory case. In the 
quantum theory, we generally deal with propagation between fixed points 
in configuration space, rather than with phase space point. Therefore, in 
the free particle case, consider the change of variables from x , po to x , x/, 
where 

x/ = x + (21) 

Hence x^ is the position after evolution for starting from x for parameter 
time t. The probability then becomes 

PA = yJ d n x f d n x w(p ,x ) 9 (j™ dt f A (4(t)) - e) (22) 
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where po = M(x/ — x )/r and 

x^(t)=x + (X/ ~ Xo) t (23) 

T 

The parameter r may in fact be scaled out of the whole expression, hence 
the probability is independent of it. 

Figure 1 

The rewritten classical probability in terms of a sum over initial and 
final points x and x^ . The probability for not entering A is a sum over paths 
as in case (a). The probability for entering A includes a sum over classical 
paths in which A lies between the initial and final points, as in case (b). But, 
to agree with the phase space form of the result ftTQ), it must also include a 
sum over initial and final points for which A does not lie between them, as 
in case (c). This figure also applies to the semiclassical propagator in the 
quantum case. 
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o ■*/ 



The result now has the form of an integral over "initial" and "final" 
points, analogous to similar results in quantum theory. The result is again 
essentially a sum over classical trajectories with the trajectories now labeled 
by any pair of points x , x/ along the trajectories, and is invariant under 
shifting x or xj along those trajectories. Naively, one might have thought 
that the restriction to paths that pass through A is imposed by summing over 
all finite length classical paths which intersect A as they go from the "initial" 
point x to "final" point x/, that is, A lies between the initial and final points. 
This is also what one might naively expect in the quantum theory version. 
However, one can see from the above construction that the correct answer is 
in fact to sum over all classical paths (which can be of infinite length) passing 
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through x and x/ that intersect A at any point along the entire trajectory, 
even if A does not lie between the two points (see Fig.l). This feature is 
related to the reparametrization invariance of the system. 

The above point turns out to be quite crucial to what follows in the 
rest of this paper, so it is worth saying it in an alternative form. Loosely 
speaking, the statement is that only the entire classical path respects the 
reparametrization invariance associated with the constraint equation. A sec- 
tion of the classical path does not. This may be expressed more precisely 
in terms of the function A(x, po,x ) introduced in Eq.([L6|). This function is 
concentrated on the entire classical trajectory, and is zero when x is not on 
the trajectory. It is easy to see that it has vanishing Poisson bracket with 
the Hamiltonian H = if(po,Xo), since we have 

/oo 
dt {if,5 (n) (x-x d (i))} 
-oo 

/OO (j 
dt -5 (n) (x-x d (t)) 
-oo dt 

= (24) 

This is the precise sense in which the entire trajectory is reparametrization 
invariant, and the phase space function A may be regarded as an observable 
- a quantity which commutes with the constraint H p3| , |2"5| . By way of 
comparison, consider a second phase space function similarly defined, but on 
only a finite section of trajectory, 

£?(x,p ,x ) = j T dt 5W(x-x d (0) (25) 
Jo 

It is easily seen that 

{H, fl(x, p , x )} = -5(x - x d (r)) + <J(x - x d (0)) (26) 

Hence B "almost" commutes with H, failing only at the end points, and 
it is in this sense that a finite section of trajectory does not fully respect 
reparametrization invariance. 

A third version of the classical result is also useful. It is of interest to 
obtain an expression for the probability for intersecting an (n— l)-dimensional 
surface S. Since the result (|l~9| ) involves the parameter time spent in a finite 
volume region A it does not apply immediately. However, suppose that 
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the set of trajectories contained in the probability distribution w intersect 
the (n — 1) -dimensional surface £ only once. Then we may consider a finite 
volume region A obtained by thickening £ along the direction of the classical 
flow. If this thickening is by a small (positive) parameter time At, then the 
quantity appearing in the ^-function in (|19"D is 

JdtJ d n x 5 (n) (x - x c/ (t)) = AtJdtJ d n ~ l x n ■ d ^^~ 

x 5 (n) (x-x c/ (t)) 
= At/[S,x d (t)] (27) 

where n is the normal to E, and we suppose that the normal is chosen so that 
n ■ dx. d /dt is positive. The quantity J[£,x d (i)], in a more general context, 
is the intersection number of the curve x d (i) with the surface S, and takes 
the value for no intersections, or ±1 (depending on whether there is an 
even or odd number of intersections). In this case we have assumed that the 
trajectories intersect at most once, hence I = or 1. We then have 

6(AtI - e) = 6(1 -e) = I (28) 

(where e = At e') and the probability for intersecting E may be written 

Pe = JdtJ d n Po d n x w(po, x ) J d n ~ l x n • <^(x - x d (t)) (29) 

At each t, we may perform a change of variables from p, x to new variables 
p' = p d (t), x' = x d (t), and using Eq . (|i~4f) , we obtain the result 

pj:= jiI dt L dnp ' dn ~ v n ■ p/ w(p/ ' x/) (30) 

Finally, the integrand is now in fact independent of t, so the t integral leads 
to an overall factor. (This might be infinite but is regularized as discussed 
below). We therefore drop the t integral. 

This result is relevant for the following reason. In the heuristic "WKB 
interpretation" of quantum cosmology, one considers WKB solutions to the 
Wheeler-DeWitt equation of the form 

* = Ce iS (31) 
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It is usually asserted that this corresponds to a set of classical trajectories 
with momentum p = VS, and with a probability of intersecting a surface £ 
given in terms of the flux of the wave function across the surface [^|, |TI| [TjJ. 
As we shall show, from the decoherent histories analysis, the quantum theory 



gives a probability for crossing a surface £ proportional to Eq . (|30D with w 
replaced by the Wigner function of the quantum theory. The Wigner function 
of the WKB wave function is, approximately HlPf , 

W(p,x) = |C(x)| 2 5(p-V5) (32) 

Inserting in Eq. (p0[) , we therefore obtain, up to overall factors, the probability 
distribution, 

p E = / d n ~ l x n • VS |C(x)| 2 (33) 
We therefore have agreement with the usual heuristic analysis. 



3 The Decoherent Histories Approach to Quan- 
tum Theory 

Our aim is to analyze the quantum case using the decoherent histories ap- 
proach to quantum theory [0. We first give a very brief review of the for- 
malism. 

In non-relativistic quantum mechanics, quantum histories are represented 
by so-called class operators Cg_, which are given by time-ordered sequences 
of projection operators 

Ca = P an (t n )---Pa 1 (tl) (34) 

(or more generally, by sums of terms of this form), where a denotes the 
string of alternatives a%, az ■ • ■ a n , and P ak {t) are projection operators in the 
Heisenberg picture. The central object of interest is then the decoherence 
functional, 

D{a,gt) =Tr(C*pC^) (35) 

Intuitively, the decoherence functional is a measure of the interference be- 
tween pairs of histories a, a'. When its real part is zero for a ^ a', we say 
that the histories are consistent and probabilities 

p(a) = D(a,a) (36) 
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obeying the usual probability sum rules may be assigned to them. Typical 
physical mechanisms which produce this situation usually cause both the 
real and imaginary part of D(a,a') to vanish. This condition is usually 
called decoherence of histories, and is related to the existence of so-called 
generalized records. 

In the non-relativistic case, for histories characterized by projections onto 
configuration space, a path integral version of the decoherence functional is 
available, and can be very useful. It has the form, 

D(a,a') = Jvx J^Vyexp (~S[x{t)] - -%(*)]) P^o,Vo) (37) 

where the sum is over pairs of paths x(t), y(t) passing through the pairs of 
regions a, a'. This is equivalent to the form (p^), (|3"5|), when the histories are 
strings of projections onto ranges of positions. Eq.(^) is a useful starting 
point for the generalization to timeless theories. 

The power of the decoherent histories approach is that it readily general- 
izes to a variety of different situations in which time plays a non-trivial role 



[|, [36|, |37|. In particular, it may be generalized to the question of interest 
here, in which the system is in an energy eigenstate and we would like to 
answer questions that do no refer to time in any way. This generalization 
requires, however, specification of the inner product used to construct the 
decoherence functional, and a prescription for the construction of the class 
operators. We consider each in turn. 



4 The Induced Inner Product 

For many situations, and especially for the analogous situation in quantum 
cosmology, the Hamiltonian has a continuous spectrum so the energy eigen- 
states are not normalizable in the usual inner product, 

(*i|* 2 ) = / d n x *;(x)* 2 (x) (38) 
A way to deal with this has been developed, and goes by the name of the 



induced inner product, or Rieffel induction |30], [L7|]. Consider the eigenvalue 
equation 

H\* EX ) = E\V EX ) (39) 
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where A denotes the degeneracy. These eigenstates will typically satisfy 



(*e>\>\*e\) = S(E-E')8 X x> (40) 

from which it is clear that the inner product diverges when E = E'. The 
induced inner product on a set of eigenstates of fixed E is defined, loosely 
speaking, by discarding the 5-function 8(E — E'). That is, the induced or 
physical inner product is then defined by 

{^E\>\^Ex) phys = 8xx> (41) 

This procedure can be defined quite rigorously, and has been discussed at 
some length in Refs.[|3C, 17]. We will use it here to construct the decoherence 



functional. A simple prescription for using it in the decoherence functional 
is to regularize each propagator and energy eigenstate by using a different 
energy for each. The final answer will then involve a number of (^-functions 
in energy, as in (|40|), which are simply dropped. 

Applied to the Klein-Gordon equation in flat space, where solutions may 
be split into positive and negative frequency parts, = + + <fi~, the induced 
inner product may be expressed in terms of the Klein-Gordon inner product 
(which is negative on the negative frequency solutions): 



)kg ~ W i9 ) KG 



(42) 



This is clearly positive definite. That is, the induced inner product effectively 
changes the sign of the negative frequency solutions to make the overall 
expression positive definite. Importantly, the induced inner product does 
not require a split into positive and negative frequency solutions, and so is 
generally applicable in quantum cosmological models. 

Although the induced inner product has been around for a long time 
(in the Klein-Gordon case it was noted by Henneaux and Teitelboim ||20|| ) , 
it is only comparatively recently that it made an appearance in quantum 
cosmology. 



5 The Class Operators 

The most important part of the construction of the decoherence functional 
for our timeless model is the construction of the class operators. On the face 
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of it, a natural prescription for constructing these is the following expression: 

/oo r 
dr e' iEr / £>x(t)exp(iS[x(t)]) 
-oo J 

x 9 ^ dt / A (x(t)) - e) (43) 

This expression consists of a functional integral between the initial and final 
points in fixed parameter time r over paths constrained to enter the region A 
(enforced by the function as in the classical case), followed by a sum over 



all possible parameter times r ||16|| . In keeping with the general ideas of Dirac 
quantization, we expect that, in order that the reparametrization invariance 
of the theory is fully respected, a properly constructed class operator should 
be annihilated by the constraint, 

HC A = (44) 

It is straightforward to show that this is indeed the case in the limit that A 
becomes the entire configuration space. However, as shown in Ref . |L5 ] , this 
is not in fact the case when the region A is finite: one obtains (^-functions on 
the boundary of A on the right-hand side of Eq.(44). 



This is a serious difficulty since it means that reparametrization invariance 
is in some sense violated. It is due to the fact that, as the end-points move 
from inside A to outside, the class operator changes discontinuously. This in 
turn is related to the fact that a projection operator onto the region A does 
not commute with the constraint equation. 

Because of this difficulty, it is necessary to replace the class operator 
Ca with a modified class operator C' A which is as much as possible, defined 
by a sum over paths passing through A but satisfies the constraint equation 
everywhere [fl7ll . At present, there does not appear to be a universally agreed 
way to do this, but some suggestions and discussion of this point were given 
in Refs.0, MM- 

Fortunately, in the semiclassical approximation, it seems clear how to 
construct the modified propagator. It is 

C A ( X/ ,x ) = 6 ^dt f A (4(t)) ~ e) P(x f ,Xo) e^f^ (45) 

Here Pe tA is the usual unrestricted semiclassical propagator, so A(x/,Xo) is 
the classical action between initial and final points, and P is a prefactor. 
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The ^-function here is the same as in the (rewritten) classical case Eq.(|22|) in 
terms of "initial" and "final" points, where Xq(£) denotes that classical path 
from x to xj. (This is exactly as in the classical case depicted in Fig.l). 
Note also that 

VA • V 6 (7^ dt f A (4(t)) - e) = (46) 

as may be shown by shifting the t integration. It follows that the modi- 
fied class operator is a semiclassical solution to the constraint equation, as 
required. 

It is important that t is integrated over an infinite range in the quantity 
inside the ^-function, otherwise the modified class operator would not in 
fact satisfy the constraint. Recall that the originally defined class operator 
Eq. (|43| ) contained a similar ^-function, with a finite range of time integration, 
which one might have been tempted to use in the semiclassical approximation, 
but this class operator does not in fact satisfy the constraint. 

Hence we see that the difference between the modified and original class 
operators in the semiclassical approximation is the difference between using 
the entire classical trajectory or using finite segments of it in the ^-functions. 
We also see that these modified class operators are the correct ones to use in 
order to be consistent with the discussion of the classical case and Eq.([22|). 
There, we saw that it is appropriate to sum over classical paths intersecting 
A even if A does not lie on the segment of classical trajectory between x 
and xj. This feature therefore appears to be necessary for the particular 
type of reparametrization invariance used here. Only the entire trajectory is 
reparametrization-invariant notion. A finite section of trajectory is not. (See 
Ref.fi~2|| for a further discussion of reparametrizations in this sort of context). 



6 Decoherence and the Decoherence Func- 
tional 

In terms of the class operators C a and an "initial" state \I/(x) (an eigenstate 
of H), the decoherence functional is 

D(a,a') = J d n x f d n x d n y C Q ( X/ , x )C a /(x / , y ) ^(x )^*(y ) (47) 
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Here, a denotes the two sets of histories that either pass through A or not. 
The decoherence functional will generally not be diagonal except for special 
initial states, so a decoherence mechanism is necessary. 

The details of the decoherence mechanism are somewhat lengthy and 



may be found in Ref.||15||. Briefly, one may couple the system to a generic 
environment consisting of large number of degrees of freedom, and then trace 
out the environmental coordinates. We do not address the more physical 
question of what the environment actually is for minisuperspace models, but 
it is sufficient to observe that there are very many fields that are unobserved 
in cosmo logical models, and so may serve as an environment. The overall 
result is that, in the semi classical approximation, the decoherence functional 
Eq.fl4"T|) becomes 

D(a,a') = J d n x f d n x d n y C a (x/, x )C af /(x / , y ) 

x F( X/ ,x ,yo) tf(x )tf*(y ) (48) 
where F(xj,x ,yo) is the influence functional in the semiclassical approxi- 



mation. As shown in Ref.||15||, it typically has the form 

1 



exp [iv a T a - ^v a v b a ab j (49) 

where v = x — y, and T a and a ab are real coefficients (and depend only on 
x + y, and a a b is a non-negative matrix. The form of the influence functional 
broadly indicates that interference between different values of x (or more gen- 
erally, different semiclassical trajectories) is suppressed, as expected. (The 
matrix a a b may develop a zero eigenvalue, but this corresponds to the fact 
that interference between points along the same trajectory is not suppressed, 
i.e., that the influence functional is sensitive only to different trajectories, 
reflecting the underlying reparametrization invariance). 

Given decoherence, we may therefore compute the probability for pass- 
ing through the region A. Introducing the variable Xo = |(x + y), and 
introducing the Wigner function, 

W(p, X) = j±pJ d n v e-* v p(X + i v , X - i v ) (50) 

it is readily shown that the probability is 

Pa — [ d n Po d n X e( f°° dt / A (X d (t)) - e) W{p , X ) (51) 
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Here, X d (t) is the classical path with initial data Xo, Po, and we have defined 
the smeared Wigner function 

W{p , X ) = J d n p exp (-ifo - p - r)v afe 1 (p - P - r) 6 ) 

xW(p,X ) (52) 



Eq. (|51f) is the main result. It shows that the decoherent histories analy- 
sis very nearly produces the classical result Eq. fll9|) , except that instead of 
a classical phase space probability distribution, we get a smeared Wigner 
function. We also see that the coefficient r a from the influence functional 
produces a small modification in the classical equations of motion, and that 
there are fluctuations about the classical motion governed by . (A pos- 
sible zero eigenvalue for simply means that the exponential in Eq.(|52j) 
becomes a 5-function in one particular direction, so there are no fluctuations 
in that direction). 



7 Summary and Discussion 

The present work shows that it is possible to carry out a decoherent histo- 
ries quantization of simple quantum cosmological models. This quantization 
method reproduces, approximately, heuristic methods based on essentially 
classical ideas. These models have the key property of not possessing a time 
parameter, but this does not appear to present an insurmountable obstruc- 
tion to the application of the approach. 

Central to both the classical and quantum problems is the notion of an 
entire trajectory. At the classical level it appears to be the appropriate 
reparametrization-invariant notion for the construction of interesting proba- 
bilities. At the quantum level, the decoherent histories approach appears to 
handle the problem in a natural way, perhaps because it readily incorporates 
the notion of trajectory. 

This approach to quantizing cosmological models in this way is certainly 
only a first bite at the problem, and a list of further topics and related is- 
sues is described at some length in Ref.[15[|. Whilst it is gratifying that 



quantum cosmological models may be successfully quantized using the de- 
coherent histories approach, what has been particularly interesting for me 
is the observation that this area appears to be a point of confluence for 
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a number of different fields. As indicated at the beginning of this article, 
for example, quantum cosmology puts severe pressure on the foundations of 
standard quantum theory. There has also been an interesting interplay go- 
ing on between quantum cosmology and the question of time in standard, 



non-relativistic quantum mechanics |3J], and this is likely to prove fruitful 
for quantum cosmology. In the long run, quantum cosmology may or may 
not be able to provide useful, testable cosmological predictions. But what is 
pretty certain is that work in this area is likely to remain a source of stimulus 
for a long time to come. 
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